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ABSTRACT 

We report on the effects of cosmic rays (CRs) on the abundance of CO in H 2 
clouds under conditions typical for star-forming galaxies in the Universe. We dis¬ 
cover that this most important molecule for tracing H 2 gas is very effectively de¬ 
stroyed in ISM environments with CR energy densities Ucr ~ (50—10 3 ) xUcR^Gai, 
a range expected in numerous star-forming systems throughout the Universe. 
This density-dependent effect operates volumetrically rather than only on molec¬ 
ular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is 
facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel 
activated by CR-produced He + . The effect we uncover is strong enough to ren¬ 
der Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus 
CO-untraceable), even in ISM environments with rather modestly enhanced av¬ 
erage CR energy densities of Ucr ~ (10 — 50) x UcR.Gai- We conclude that the 
CR-induced destruction of CO in molecular clouds, unhindered by dust absorp¬ 
tion, is perhaps the single most important factor controlling the CO-visibility of 



2 


molecular gas in vigorously star-forming galaxies. We anticipate that a second 
order effect of this CO destruction mechanism will be to make the H 2 distribution 
in the gas-rich disks of such galaxies appear much clumpier in CO J=1-0, 2-1 
line emission than it actually is. Finally we give an analytical approximation of 
the CO/H 2 abundance ratio as a function of gas density and CR energy density 
for use in galaxy-size or cosmological hydrodynamical simulations, and propose 
some key observational tests. 

Subject headings: (ISM:) cosmic rays - astro chemistry - methods: numerical - 
ISM: abundances - (ISM:) photon-dominated region (PDR) - galaxies: ISM 


Introduction 


Using the mm/submm rotational transitions of CO, the most abundant molecule after 
Ho (CO/H 2 ~ 10 -4 ), as tracers of the Ho distribution and average physical conditions in 
galaxies is now a well-established method, and indeed a main driver behind the largest 
ground-based astronomical facility on Earth, the Atacama Large Millimeter Array (ALMA) 
in Llano Chajnantor of North Chile. Ever since the first detection of a bright CO J = 1 — 0 


line in Orion (Wilson et ah 1970) there has been tremendous progress towards detecting 


ever higher-J transitions of CO, in our own Galaxy as well as in others (e.g. Cronin et al. 


1976 Guesten et al. 1993 Devereux et al. 1994; Dumke et al. 

2001; Yao et al. 

2003a ; b Mao 

et al. 2010 

Papadopoulos et al. 2007 Papadopoulos 2010 

Leech et al. 2010 

) and of other 

molecules such as HCN, HCO + , and CS (Nguyen et al. 1992| 

Solomon et al. 

1992f Paglione 

et ah 1997 

Gao & Solomon 2004 Gracia-Carpio et al. 2006, 

2008 Papadopoulos||2007; Krips 


2008 

Greve et al. 

2009; 

Zhang et al. 

2014) 


then be used as powerful probes of the average density, temperature, and dynamical state 
(i.e. gravitationally bound or unbound) of the molecular gas in galaxies (Papadopoulos et al. 


2012, and references therein). A lot of effort has been focused particularly towards using the 


two lowest CO rotational lines as H 2 gas mass tracers in galaxies via the so-called A'ccrfactor 


(e-g- 

Dickman et al. 

1986 

Maloney & Black 

1988 

Young & Scovillc 

1991 

Solomon et al. 

1992 

Bryant & Scovillc|l996 

Solomon et al. 

L997; 

Downes & Solomon 1998) 

In this context 


higher-J CO and other molecular lines remain very important since the constrains they set 
on the average density, thermal and dynamical state of H 2 cloud ensemble^ help determine 
A"co and find whether it deviates from its Galactic calibrated value A"co,Gai ~ 5 A-fl (see 


1 A'co retains its mass-tracing value only for cloud ensembles rather than cloud parts or single H 2 clouds 

2 X( = M 0 (K kms -1 pc 1 2 ) -1 
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Papadopoulos et al. 2012, for a recent such study) 


Nevertheless, all available studies quantifying the effects of the H 2 physical conditions 
on the (H 2 mass)/(CO luminosity) ratio rely on the premise of a CO-rich H 2 gas phase. 
Significantly lower CO/H 2 ratios would render the use of CO lines as molecular gas mass 
tracers problematic. Extensive observational and theoretical work has indeed shown that 
a CO-poor phase becomes possible for ISM with low average metallicities (Z< 0.2 Z Q ) and 


strong FUV fields (e.g. Maloney & Black 1988 Madden et al. 1997 Israel 1997 Pak et al. 


1998 Bell et al. 2006a). Metallicity (Z) and the average FUV radiation held (G 0 ) then 


emerge as the two most important factors determining CO/H 2 in metal-poor ISM, while 
models that incorporate the underlying ISM physics have been made to correct for such 


CO-dark H 2 mass per cloud (Pak et al] 1998 Bolatto et al. 1999). These can be of help as 


long as: a) some CO survives in the inner cloud regions (so that CO line emission can still 
mark such clouds), b) the average metallicity and mean FUV radiation held are known (so 
that their values can be inserted in the models). Such models will obviously not work for 
metallicities that fall too low and/or average FUV helds that are strong enough to keep CO 
fully dissociatec0 (e.g. Bolatto et al. 1999 Bell et al. 2006b). 


Despite existing theoretical work that investigated the effects of CRs on ISM chemistry 


e.g. Bayet et al. 2011 Meijerink et al. 2011), the effects on CO/H 2 of CRs (which, unlike 


FUV radiation, travel nearly unimpeded through the irrespective of its dust content), or 
their impact on the observability of H 2 gas via CO lines has not been investigated in detail. 


Recently, Bialy & Sternberg (2014) presented an analytical study of the importance of CRs in 


the abundance of the most important species under different metallicities, using isothermal 
calculations. In this work, we present our investigation of the effects produced by CRs 
starting with the ISM of metal-rich star-forming galaxies but, unlike Bialy & Sternberg 
(2014), we perform full thermochemical calculations (i.e. we solve for thermal balance) 


while using a much more extensive chemical network. We assume CRs fully penetrating the 
molecular clouds, with an average CR energy density scaling as Ucroc Psfr (where Psfr 
is the average star-formation rate energy density; see |Papa dopoulos 2010, and references 
therein for detailed arguments). 


3 For clarity we must emphasize here that even in the CO-rich parts of the Galaxy CO-poor regions do 
exist (e.g. diffuse H 2 clouds, outer parts of GMCs, regions very near O, B stars), but their masses are 
typically too small to make any impact on the total H 2 mass budget as traced by the CO-rich H 2 gas. 
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2. The effect of CRs on the CO/H 2 ratio in molecular clouds 


We first perform a grid of runs of one-dimensional uniform density clouds in which we 
vary the density, Ur, from (10 2 — 10 5 ) cm -3 and the cosmic ray ionization rate, Ccr ; from 
10 -1 ' — 10 -13 s -1 . The choice for the Ccr range is dictated by that deduced within the Milky 
Way (~ (1 — 5) x 10 -17 s -1 ), and from the Ccr expected in galaxies with star-formation (SF) 
rate densities of psfr~ (1 — 10 3 ) x psFR.Gai (assuming Ccr oc Ucr where Ucr is the CR 
energy density, and Ucr oc Psfr for a CR-transparent ISM). 


In order to isolate the effects of CRs on the important CO/H 2 ratio we set the incident 
FUV held G = 0 G a (G 0 the average interstellar radiation held as described by Drain e|l9 78), 
and consider only Solar metallicity (Z=Z Q ) gas. We chose a density range that encompasses 
the low-density Galactic regions where the Hi—>-H 2 phase transition takes place in the Cold 
Neutral Medium ((uh)~200 cnr 3 ), to typical Galactic GMCs ((nn)~500 — 10 3 cm -3 ), and 
up to the high densities typical of SF cores ((uh)~ 10 5 cm -3 ). We then use the 3d-pdr code 
(see Appendix [A]) to conduct full thermochemical calculations on a uniform density cloud 
permeated by an average Ucr (that proportionally sets Ccr) for the grid of [ur — Ccr] values. 
The results on the CO/H 2 ratio (blue lines, marking CO/H 2 = 10 -5 ) can be viewed in Fig. 
[lj where we also show the Cl/CO ratio (green lines, indicating the Cl/CO=5 ratio) since 
atomic carbon is the dominant species resulting from the CR-induced destruction of CO. 


From Fig|l] it can be seen that molecular clouds with densities typical of the Galactic 
GMC density range (~10 2 — 10 3 cm -3 ) can become very CO-poor even in ISM environments 
with SFR densities (and thus Ucr) only ~ 10 x Galactic. This is even more the case for the 
ISM of gas-rich SF systems like the so-called BzK galaxies. These are isolated disks, typically 
at z~ 1 — 2, with SFR ~ (50 — 150) x Galactic, and representative of a major SF mode in 
the Early Universe (Daddi et al. 2010). In such systems CO/H 2 < 10 -5 (10 times less 
abundant than in the Galaxy) for much of the range of average H 2 gas densities typical for 
Galactic GMCs. Thus, unless the molecular gas in such galaxies lies preferentially in cloud 
structures with n(H 2 ) > 10 3 cm -3 , CO imaging for the purpose of revealing H 2 gas disk sizes, 
average molecular gas physical conditions, and H 2 mass distribution (e.g. Daddi et al. 2010 


2014) may be inadequate. Figure [I] also reveals that efficient CR-induced CO destruction 
encompasses ever increasing gas densities for higher Ccr> reaching up to n 10 4 cm 3 for Ccr 
values expected in compact merger/starburst galaxies in the local Universe (ULIRGs), and 
similar systems like the so-called Submm Galaxies (SMGs) found at much larger numbers in 
the distant Universe (Smail et al. 1997 Hughes et al. 1998) 


Furthermore, even in the Galaxy there seems to be a [?zr — Ccr] parameter space where 
low-density molecular gas (n H ~ (1 — 3) x 10 2 cm -3 ) can be very CO-poor. Such a gas 
phase may indeed exist in the outer Galactic disk (i?cai > 15kpc) where a combination of 
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Fig. 1.— Map showing the traceability of H 2 using CO and Cl as a function of H number 
density % and the cosmic-ray ionization rate (cr- The three horizontal areas correspond 
to the £cr range for (from bottom to top) the Milky Way, BzK galaxies, and Starburst 
galaxies. The shaded area below the blue line (CO/H 2 = ICC 5 ) corresponds to the region for 
which CO and Cl can both trace H 2 . The region above the blue line and below the top left 
dark-shaded region corresponds to the area in which CO has been destroyed by cosmic-rays 
but H 2 survives, making Cl the only possible tracer. The dark-shaded region on top left 
corresponds to the area where H 2 is also destroyed by the interaction of cosmic rays hence 
being rich in Hi (Hi/ 2H 2 = 1). Note that for £cr > 10 -14 s -1 CO is destroyed less effectively 
for lower densities than % = 10 4 cm -3 due to the destruction of He + . In 14.1 we discuss the 
chemistry behind this behaviour. 
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low metallicities (Z ~ 0.15 — 0.2 Z & ) and low densities keep CO formation rates low while a 
still substantial FUV radiation field maintains high CO destruction rates. This can yield a 
CO-dark H 2 phase that can be quite massive (Papadopoulos et al. 2002). Any additional CO 
destruction due to CRs (which propagate freely to the outer Galactic disk) can only further 
reduce CO/H 2 , rendering the corresponding molecular gas phase even harder to detect via 
CO line emission (see also Wolhre et al. 2010, for a recent study of CO-dark H 2 gas in the 
Galaxy). 


2.1. The intensity of CO J=l—0, 2—1 and Cl 1—0, 2—1 lines 


We will now examine the emergent intensities of the CO J=l-0, 2-1 (~115GHz, 
~230 GHz) and Cl J=l-0, 2-1 (~492 GHz, ~809GHz) lines across the [nn — Ccr] parameter 
space shown in Fig. [lj in order to access their detectability in CO-poor versus CO-rich gas 
phases. We choose only the two lowest- J CO lines as the only CO transitions that are not 
excitation-biased, when it comes to the underlying state of the gas necessary to significantly 
excite them. Higher-J CO lines will trace mainly warm and dense H 2 gas (n > 10 4 cm -3 , 
Ej/ke > 30 K), and are not reliable tracers of the general H 2 gas distribution irrespective of 
its state. Thus J=l-0, 2-1 are the CO lines where most observational calibrations of the so- 
called X co =M(H 2 )/L co factor have been performed (see Bolatto et al. 2013, for a review). 
The two Cl lines on the other hand are powerful alternatives for tracing the general H 2 gas 
distribution (especially Cl J=l-0), often being brighter than CO J=l-0, 2-1 in CO-rich gas 
while remaining luminous even in CO-poor gas (Papadopoulos et al. 2004). 


Figure [2] plots the line-integrated intensities of the CO J=l-0 and J=2-l transitions. It 
can immediately be seen that these intensities fall precipitously towards the CO-poor area. 
Indeed, their strengths relative to those typical in the CO-rich regions range from ~0.39 at 
the “crossing” of the blue line (used to mark the CO-rich/CO-poor transition zone in the 
[n H — Ccr] parameter space), to only ~0.025-0.15 deeper into the CO-poor regions in the 
upper left corner of the parameter space. Thus CO line observations conducted to yield 
the same S/N ratio for CO-poor molecular gas distributions would need ~7-1600 longer 
integration time with respect to those of CO-rich regions. The Cl J=l-0, 2-1 lines on the 
other hand shown in Fig. [3j remain bright throughout most of the [nn — Ccr] parameter 
space, and actually much brighter than the two CO linc^] This enables the two Cl lines to 


4 This does not translate to Ici/Ico x higher S/N ratios for Cl versus CO observations as T sys (v) will be 


higher at Cl line frequencies. Nevertheless, significant gains in S/N are expected (e.g. Papadopoulos et al. 


2004). 
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trace H 2 gas both in CO-poor and CO-rich regions. 

There is also an area of low CO and Cl line intensities in the (high-nn)/ (Iow-Ccr) region 
of the [iih — Ccr] parameter space (lower right corner) in all plots of both Figs. [2] and [3} This 
is the area where non-SF, dense gas cores embedded deep into Galactic GMCs are located 
(e.g. Bok globules). These are very cold regions (T gas ~10K), and the reduction of CO 
<7=1—0, 2-1 line intensities with respect to the warmer gas in the (lower-n H )/(higher- Ccr) 
areas almost purely reflects a temperature drop (since CO <7=1-0, 2-1 are thermalized and 
optically thick with Ei 0 , 2 i/(kBT gas ) < 1). The drop of the Cl line intensity towards that 
area is due both to this temperature drop (that contributes non-linearly to the intensity 
reduction since for these lines E 1021 /(k B T gas ) > 1, especially for Cl 7=2-1), and a [C/CO] 
abundance reduction. Nevertheless, given that such dense, cold gas cores contain only small 
mass fractions (<1%) of Galactic GMCs, we do not anticipate that this (CO and Cl)-dark 
region of [n B — Ccr] parameter space to be of any consequence when it comes to the inventory 
of H 2 gas mass. 


CO 1-0 



CO 2-1 



Fig. 2.— Logaritmic contour /grayscale plots showing the brightnesses of CO 7=1-0 (left) 
and CO 7=2-1 (right) in the [n H — Ccr] space. The brightness is calculated by integrating 
the local emissivity for a depth of Ay = 10 mag, and is in units of erg cm -2 s -1 sr . The 
blue line corresponds to CO/H 2 = 10 -5 (see Fig. [l] for further details). The large reduction 
of the CO line strengths in the CO-poor regions is evident (see §|2.1|for detailed discussion). 


2.2. The effect of higher FUV fields 

The choice of G = 0G o in our calculations was driven by our intent to isolate the effect 
of CRs on the average [CO/H 2 ] abundance in molecular clouds. However ISM environments 
with higher SFR densities boost both the average Ccr an d the average FUV field G. Nev- 





















Cl 1-0 



Cl 2-1 



Fig. 3.— Logaritmic contour/grayscale plots showing Cl J=l-0 (left) and Cl J=2-l (right) 
brightness in the [n H — Ccr] space. The brightness is calculated by integrating the local 
emissivity for a depth of Ay = 10 mag, and is in units of ergcm _2 s _1 sr _1 . The blue line 
corresponds to CO/ff 2 = 10~ 5 (see Fig. [l] for further details). The two Cl lines remain 
luminous for much of the [nn — Ccr] space, even in CO poor regions. They become faint (as 
well as the CO ,7=1-0, 2-1 lines) only for dense/cold gas with Galactic levels of Ccr (see 


12.1 for detailed discussion). 


erthcless, unlike Ccr the boost of the average FUV radiation field is much reduced by the 
dust absorption in the metal-rich ISM environments we study here. Indeed it is expected to 
be caped to G ~ 10 3 G o (Papadopoulos et al. 2014), with higher values such as ~ 10 4 — 10 5 


associated with gas regions near O, B stars that contain only small mass fractions of molec¬ 
ular clouds. Furthermore, in environments of high gas turbulence (typical for ISM with high 


SFR densities, see also discussion in 15.1) such clouds are expected to have large ( Ay ), much 
larger than (Ay) ~ 7 — 8 found for Galactic GMCs (see Pclupessy et al. 2006). Thus we do 


not expect the [CO/H 2 ] to be affected by the average FUV fields for much of the molecular 
gas mass, a state of affairs which may, however, change for low-metallicity ISM. We must 
however note that increasing the FUV will result to a higher Cll abundance (and Cll line 
emission) in low densities. 

We nevertheless still performed calculations for G = 10 3 G o and G = 10 5 G o , and found 
that the corresponding map shown in Fig. [I] to be in excellent agreement for Ay > 7 mag. 
In addition, we compare the gas temperature of selected simulations with those studied in 


Meijerink et al. (2011) who included the effect of FUV radiation and also find excellent 


agreement. 
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2.3. 


Comparison with other theoretical models, and a [CO/H 2 ] sensitivity on 

T 

-L O 


-gas 


As we noted in the Introduction, theoretical models by Bayet et al. (2011); Meijerink 


et al. (2011); Bialy & Sternberg (2014) have examined the effects of CRs in uniform density 


clouds, even if they did not focus on the CR-controlled [CO/H 2 ] abundance and its obser¬ 
vational consequences in various extragalactic environments. Figure [4] shows a comparison 


between Bayet et al. (2011); Bialy & Sternberg (2014) and 3d-pdr results. Red and green 


lines compare 3D-PDR against Bialy & Sternberg (2014) respectively for isothermal runs at 
T gas = 100 K interacting with £cr = 10 - 16 s -1 . It can be seen that the agreement is very 


good. We note that Bialy & Sternberg (2014) and the present work used different chemical 
network which can explain the discrepancy of the two curves for tt-h < 2 x 10 3 cm -3 , while for 
«h > 2x 10 3 cm -3 the agreement is excellent. Switching on the thermal balance calculations 
(blue line) makes the CO/H 2 ratio lower for % < 2 x 10 3 cm -3 while for ~ 10 2 cm -3 this 
ratio is < 10~ 5 , differing by an order of magnitude with respect to the isothermal runs. This 
demonstrates the importance of thermal balance calculations in estimating the effect of CRs 


on CO/H 2 as the latter appears to have high sensitivity on gas temperature (see also 12.6). 
Lower temperatures form 0 2 via the reaction 


O + OH -> CL 


H 


( 1 ) 


which then reacts with C to form CO as shown in Reaction 12 (see 14.1), hence increasing 
the abundance of the latter. 


The above reaction rate can be estimated using the following equation (e.g. Millar et 
al. |1997): 




( 2 ) 


with a = 3.69 x ICC 11 , /3 = —0.27 and 7 = 12.9. For T gas = 14 K (as obtained by the 3d-pdr 
runs), k = 3.36 x 10~ n while for T gas = 100 K, k = 4.36 x 10 -11 . Hence at lower temperatures, 
such as those found by the thermal balance iterations in this work, the abundance of CO 
will be higher. 


Finally, the yellow circular point corresponds to a simulation presented in Bayet et al. 


( 2011 ) who included thermal balance calculations but used different chemical network and 


initial elemental abundances of species resulting in a discrepancy between their work and 


the present paper. By using the Bayet et al. (2011) initial elemental abundances, we obtain 
excellent agreement. 
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n x [cm ! ] 


Fig. 4. — Comparison between 3d-pdr, 
The cosmic-ray ionization rate is £cr = 10 



& Sternberg (2014) isothermal simulations and the red line mark our 3D-PDR isothermal 
calculations at T gas = 100 K. The blue line corresponds to the 3D-PDR calculations that 
solve also for gas thermal balance. The yellow circular point corresponds to a model by 

See 


Bayet et al. (2011) for % = 10 4 cm 


3 


for the relevant discussion. 


2.4. A low-density CO/H 2 upwards turnover at high Ucr values 


The CO/H 2 = ICR 5 contour marks a near-linear path across the [nn — Ccr] parame¬ 
ter space for up to Ccr ~ 10 _14 s _1 . However, for higher Ccr values there seems to be a 
“turnover” of this CO/H 2 contour value towards lower densities, i.e. low-density gas starts 
becoming CO-rich again despite the now exceptionally high CR energy densities. We give 


a full description of the underlying chemistry of this behaviour in §4.1| We must note that 
although such CR-intense environments are too extreme to be considered representatives 
of average conditions in the ISM of SF galaxies, they could be possible locally near highly 
energetic CR-accelerating phenomena such as supernova explosions and young SNRs. 


2.5. A non-uniform cloud model embedded in a 1G 0 FUV field 

While our investigation using uniform density 3D-PDR models allows a clear overview of 
the density-dependence of the CR-induced destruction of CO, and what type of extragalactic 
environments would be most affected (see Fig |T|) , it is obviously too simple of an approach 
when it comes to individual molecular clouds where hierarchical density structures are known 
to exist. For the present work we choose to investigate the CR-controlled CO/H 2 and 
Cl/CO ratios for a cloud with a simple nn(r) oc r ^ 1 profile. In the future we will use 
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the full capabilities of our 3d-pdr code where fractal density structures have been recently 
incorporated to examine the CR effects on the CO and Cl abundances in much more realistic 
settings (i.e. such as those examined by Walch et ah 2012). 


We consider a cloud with radius R = 20 pc embedded in an environment with different 
cosmic-ray ionization rates. In particular, we perform a suite of runs in which we use a 
density profile following the Larson ( |1981 ) power law of nu(r) oc r~ 1 by adopting the density 
function 


n H (r) 


4 x 10 3 


0.04 + 



cm 


(3) 


The above density distribution has the following properties: the central core of the cloud 
reaches densities up to tt-h = 10 5 cm -3 while at R = 20 pc it reaches densities down to 
Tin = 200 cm -3 . The total mass of the cloud is Mtot = 2.47 x 1O 5 M 0 while its average 
density is {%) — 300 cnr 3 , typical of Galactic GMCs. We also include an external FUV 
radiation held of intensity G = 1G 0 . Furthermore, the visual extinction at the centre of the 
core (measured from the edge of the sphere) is Ay, max ~ 50 mags. 

Figure [5] shows the results of the 3D-PDR calculations for different cosmic-ray ionization 
rates (£cr = 10 -17 — 10^ 14 s' 1 ). We can now see the effects of enhanced CRs within an 
inhomogeneous H 2 cloud. The red line marks a Cll/Cl=5 ratio, with the spherical shell 
outside that line (i.e. outer part) consisting almost entirely of ClI. The green line corresponds 
to Cl/CO=5 hence the shell between the red and green lines is Cl-dominated while the 
rest of the cloud (i.e. inner to the green line) is CO-dominated. The blue line defines the 
CO/H 2 = 10~ 5 ratio and hence the maximum volume where CO can easily trace H 2 . Finally, 
we also plot the Hl/2H 2 =1 ratio with a black line (visible only in Panel D). The volume inner 
to that line is H 2 -dominated, with the outer part where CRs have mostly destroyed H 2 to 
form Hi. 


In Table[j]we list the H 2 gas mass fractions traceable by CO, Cl and Cll in each case. We 
find that Cll cannot trace more than 1% of the cloud in any of the cases. As Ccr increases, the 
Ho gas becomes Cl-rich as CO gets destroyed by the reaction with He + . At £cr = 10^ 16 s^ 1 
only ~ 33% remains CO-rich while for £cr > 10~ 15 s^ 1 most of the cloud is CO-poor (but 
Cl-rich). Evidently, cosmic rays play a major role in controlling the CO-visibility of H 2 gas 
while leaving the Cl lines as the only alternative molecular gas mass tracers. 
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Fig. 5.— The spherically symmetric one-dimensional simulation for the cloud with density 
profile as described in Eqnj3j Each different panel corresponds to different £cr- From A to D, 
Ccr = 10” 17 , 10 -16 , 10~ 15 , 10 -14 s -1 . The colour bar corresponds to the «h number density. 
The reel line corresponds to Cn/Ci=5. The blue line corresponds to CO/H 2 = 10~ 5 , whereas 
the green one to Cl/CO= 5. The black line (which is only visible in the D panel because in 
all other cases it lies in the edge of the cloud), corresponds to Hl/2H 2 = 1. H 2 is traceable 
by CO and Cl in the region below the blue line; by Cl only in the region between the red 
and the blue line; and by ClI only in the region outside the red line. For the particular case 
of panel D, H 2 is traceable by Cl only in the inner area of the cloud. 
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2.6. The importance of thermal balance 


We will now briefly explore the effects of thermal balance on the results of this work. To 
do this, we switch off the thermal balance iteration scheme in 3D-PDR and instead we run 
two isothermal grids of runs. Figure [6] plots the CO/H 2 ratio in the n H and Ccr parameter 
space for T gas = 20K (left panel) and T gas = 100K (right panel). The thick solid line 
corresponds to CO/H 2 = 10 -5 . In both cases, the CO/H 2 ratio follows a linear behaviour 
with log ?t.h and log£cR thus demonstrating that the dependence of the chemical reactions 
on the gas temperature in such cosmic-ray dominated environments is a highly non-linear 


and at the same time a very important one. Indeed Bialy & Sternberg (2014), which used an 


analytical approach (and simpler chemical pathways) along with an isothermal assumption 
(at T gas = 100 K), do not recover this behavior, nor do they find the large reduction of 
CO/H 2 at high CR energy densities. 




Fig. 6.— The CO/H 2 ratio for the isothermal runs discussed in 12.6 The left panel shows 


the resultant contour plot for T gas = 20 K and the right panel for T gas = 100 K (c.f. Bialy & 
Sternberg 2014). The thick solid line corresponds to CO/H 2 = 10~ 5 . A strong dependency 


on the gas temperature is apparent, with the highest temperatures corresponding to much 
smaller surpressions of the CO/H 2 abundance ratio. This demonstrates the importance of 
solving for the detailed thermal balance. 


3. A secondary effect: clumpy low-J CO emission in gas-rich SF disks? 

A secondary effect that may occur in H 2 -rich galactic disks permeated by high Ucr 
backgrounds (boosted by the intense star formation typically hosted by such disks) which 
destroy CO more effectively at lower gas densities is to make such disks appear clumpy in 
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low-J CO line emissioif\ This can occur simply because the CR-induced destruction of CO, 
being more effective at lower densities, will leave only smaller and denser parts of molecular 
clouds to be CO-rich. The appearance of such clouds or cloud ensembles via their low-J CO 
line emission will then be much clumpier, depending on the average Ucr background and 
the underlying mass-density dM gas (n)/dn function of the molecular clouds. Finally, if lower 
density molecular clouds lie preferably at larger galactocentric distances, such CR-inundated 
molecular gas disks could appear more compact in low-J CO emission than they actually are. 
Such large scale CR-induced effects on the (CO-line)-visibility of the H 2 gas distribution in 
galaxies can only be accentuated once metallicity gradients (with lower-metallicity gas found 


preferably at larger galactocentric distances) are also taken into account (see also £5.1) 


There are now several deep CO J = 1 — 0, 2 — 1 images available for gas-rich SF 


Daddi et al. 

2010 ; 

Ivison et al. 

2010 

Riechers et al. 


2011 Hodge et al. 2012), with clumpy gas disks often claimed (e.g. Hodge et al. 2012). In 
such cases it would be interesting to also conduct deep Cl J=l-0, 2-1 line imaging since Cl 
is a by-product of the CR-induced CO destruction which remains abundant also in CO-rich 
molecular gas. Thus it can then trace any CO-dark H 2 gas along with its CO-rich phase (see 
FigJTJ) , and may thus reveal smoother H 2 gas distributions in CR-permeated SF disks than 
the CO J=l-0, 2-1 lines. 

On the theoretical front, hydrodynamical simulations of gas-rich disks with: a) high 
enough resolution to track the molecular cloud turbulence (which sets their dM gas (n)/dn) 
and b) detailed H 2 , CO formation/destruction chemical networks, can elucidate whether 
CRs can induce such a bias on the structural characteristics of H 2 -rich disks derived from 


CO J=l—0, 2-1 lines (see also £5.2). 


4. A chemical and analytical determination of the CO-dark H 2 domain 

A main result of this study is the discovery of a large area in the [rin — Ccr] parameter 
space where gas is mostly molecular but CO-poor. In this Section we identify the responsible 
chemical pathways for this behaviour, then attempt to find an analytical equation that could 
be used in purely (hydro-)dynamical codes for determining the average CO/H 2 ratio as a 
function of gas density and Ccr in turbulent molecular gas clouds. 


5 Molecular gas disks in CO J = 3—2 and higher-J CO line emission are expected to appear much clumpier 
than in CO J = 1 — 0, 2 — 1 because high-J CO lines sample only the more compact dense and warm gas 
regions near SF sites while CO J=l-0, 2-1 are expected to trace the entire molecular gas distribution (as 
long as it is CO-rich). 
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4.1. The chemistry of H 2 and CO 

While the main chemical formation and destruction routes of both H 2 and CO in the 
gas phase are very well known, in order to understand the region in the [-n H — Ccr] plane of 
Figjl] where CO is destroyed but H 2 is not, one needs to take into consideration several more 
chemical reactions as a function of these two parameters. We perform a chemical analysis 
of three representative such models outlining that region (their properties in tabulated in 
Table [2] and the rates used for the reactions presented below are shown in Table [3]). We 
note, however, that the chemical pathway and particularly the efficiency of the reaction 
rates follows a non-linear behavior, therefore the analysis below is very model-dependent 
and hence should be taken as qualitative. 

In Ml and M2, H 2 is mainly formed by the ion-neutral reaction: 

H + H+ ->■ H 2 + H + (4) 

which contributes ~(70-90) % to the formation of H 2 . For the formation of H 2 the behaviour 
is similar for the model M3, with the difference that the ion-neutral reactions forming H 2 
are: 


H+ + 0 ->■ OH+ + H 2 

(5) 

H+ + C -> CH+ + H 2 

(6) 

Hg + CO ->• HCO+ + H 2 

(7) 


Ho is then destroyed in almost equal efficiencies by recombination and cosmic rays ionization 
in all models: 


H 2 + c.r. ->• H^" + e (8) 

H+ + H 2 -> H+ + H (9) 

Hence cosmic rays, although an efficient mean of destruction, are not necessarily the dom¬ 
inant route. CO on the other hand is formed mainly via neutral-neutral reactions in Ml, 
with the two most efficient (~ 30% each) routes being: 

CH + O —)• CO + H (10) 

C + OH —>■ CO + H (11) 

with some contribution (~ 15%) from: 


C + 0 2 ->■ CO + O 


(12) 


16 - 


Table 1: This table shows the amount (%) of H 2 gas which is traceable by CO, Cl and ClI 


for the four simulations presented in 12.5 The first column corresponds to the panels of 


Fig|5j The second column corresponds to the intensity of CRs and the third column to the 


H 2 gas mass. The rest of columns correspond to the amount of H 2 that is CO, Cl, and Cll 
bright i.e. traceable by these elements. It can be seen that once we increase the Ccr value 
Cl becomes the most important tracer for H 2 . 

Panel Ccr [xHT 1 ^" 1 ] H 2 [xlO 5 M 0 ] CO & Cl (%) Cl (%) Cll (%) 


A 

1 

2.47 

~55 

~44 

<1 

B 

10 

2.47 

~33 

~76 

<1 

C 

10 2 

2.47 

<1 

~99 

<1 

D 

10 3 

1.38 

<0.1 

>99.9 

N/A 


Table 2: Properties of three representative models marking the outlines of CO-poor/H 2 -rich 
region. 


Model ID 

log (cm 3 ) 

log Ccr (s -1 ) 

Ml 

2.5 

-16 

M2 

2.5 

-15 

M3 

3.5 

-15 
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Table 3: Reaction rates for the chemical analysis discussed in £4.1 


Reaction a /3 7 


4 

6.40 

X 

10 - 

10 

0.0 

0.0 

5 

7.98 

X 

10 - 

10 

-0.16 

1.4 

6 

2.00 

X 

10 

-9 

0.0 

0.0 

7 

1.36 

X 

10 ' 

-9 

0.0 

0.0 

8 

1.20 

X 

10 - 

17 

0.0 

0.0 

9 

2.08 

X 

10 

-9 

0.0 

0.0 

10 

6.02 

X 

10 - 

11 

0.1 

-4.5 

11 

1.00 

X 

10 - 

10 

0.0 

0.0 

12 

5.56 

X 

10 - 

11 

0.41 

-26.9 

13 

7.50 

X 

10 - 

10 

0.0 

0.0 

14 

1.60 

X 

10 ' 

-9 

0.0 

0.0 


In M2, the Reaction 12 is not important but the rest of CO is formed via the ion-neutral 
reaction: 

H + CO + ->■ CO + H + (13) 

Interestingly, in M3 CO forms only via the neutral-neutral reactions listed above (i.e Reac¬ 
tions 


10 12 ). 


However, in all models CO is always very efficiently destroyed indirectly by cosmic rays, 
via He + (which forms via cosmic ray ionization) with the reaction: 


He H 


CO -)■ O + C H 


He 


(14) 


It therefore seems that a main reason for a CO-dark H 2 phase is that a higher cosmic ray 
ionization rate drives a higher rate of destruction of CO, but has a much less dramatic effect 
on H 2 . It is worth noting here that He + can be destroyed by reacting with H 2 ; this reaction 
has a small barrier leading to its rate coefficient becoming a factor of two higher at 100 K 
(i.e. the gas temperature that Bialy & Sternberg 2014;, considered) than at 40 K (i.e. the 
gas temperature of the M3 model), possibly reducing He + enough to allow CO to survive at 
higher abundances. A thorough investigation of the chemistry, as a function of temperature, 
that can lead to a CO-poor/H 2 -rich environment is beyond the scope of this paper but it is 
clear that chemistry resulting from a high cosmic ray ionization rate when thermal balance 
is taken into account is not easily determined a priori. 

Fig |l] also shows a high £cr regime with a small parameter space where the general trend 
of (higher-CcR,)— ^(higher -11 gas becoming CO-poor) is reversed. In the left panel of Fig|7]we 
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now show how CO/H 2 varies with increasing Ccr for three different number densities, while 
the right panel shows the corresponding gas temperatures. The left panel of Figj7[ for a 
density range % = 10 3 — lCrcm -3 and £cr > 3 x 10 _15 s , shows the reversal of the 

aforementioned trend, with high-density gas becoming CO-rich again despite the now very 
high Ccr values. This reversal is much less prominent for low-density gas. We ran a chemical 
analysis for these models and we found that He + is efficiently (by > 50%) destroyed by: 

He + + M —)• M + + He (15) 


where M stands for metals for which we used a total elemental fractional abundance of 
4xl0~ 5 (jAsplund et al. 2009). Note that Reaction 15 is not very efficient in models where 
CO decreases with cosmic ray ionization rate. 




Fig. 7.— The left panel shows the abundance of CO (relative to hydrogen) versus the 
cosmic-ray ionization rate (Ccr) f° r three different densities (% = 10 3 cm -3 solid line, tt-h = 
10 3 ' 5 cm -3 dashed line, and tt-h = 10 4 cm -3 dot-dashed line). For £cr > 10 _ 14 s _1 we observe 
an increase of CO abundance due to the destruction of He + (see 14.1). The right panel shows 
the corresponding gas temperature which increases with increasing ()cr. 


Nevertheless such extreme CR-permeated gas regions where this effect can take place are 
unlikely to be representative of average ISM conditions even in the most extreme merger/starbursts 
known, yet they could be found near or within young SNRs. In such a case the effect de¬ 
scribed here may indicate the survival of the CO molecule in dense molecular gas within 
some truly extreme CR-dominated regions. 
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4.2. An analytical approximation of the the CO-dark H 2 domain 

The left panel of Figj8]is a contour plot showing 3 d-pdr results of the log(CO/H 2 ) as 
a function of % and £cr- The following function is best fit to our data: 

log(CO/H 2 ) = —110.7 + 28.69 log uh 
— 12.82 log £ C r — 1.512(log n H ) 2 
+2.642 log n H logCcR - 0.5377(logC CR ) 2 
—0.08047(log n H ) 2 log Ccr 

+0.06209 log n H (log Ccr ) 2 - 0.008312(logCcR) 3 (16) 

and is plotted in the right panel of Fig|8j This equation is applicable for negligible FUV 
radiation within a cloud, typically for Ay > 7 mag (the case for most H 2 gas mass in star- 
bursts). In cases where FUV radiation is significant this parametric equation still provides 
a maximum CO/H 2 abundance ratio possible for a given average SFR density (since FUV 
can only destroy CO further while leaving the mostly self-shielding H 2 intact). 




Fig. 8.— The left panel of the figure shows a contour plot of the 3D-PDR results for the 
CO/H 2 ratio. The thick solid line corresponds to CO/H 2 = 10~ 5 . The right panel shows the 
corresponding best fit function described in Eqnfl6| 

This relation cannot track the behaviour of CO/H 2 as extracted from the detailed 3 d- 
PDR models in some domains (e.g. at % ~ 1000cm -3 and very high cosmic ray ionization). 
However for the [nn — Ccr] parameter space where most of the molecular gas mass of SF 
galaxies is expected to be this relation gives reasonably good result^] It can then be used in 


6 The coefficients have a 95% confidence bounds, and a Mean Squared Error (MSE) analysis on this 
function yields an MSE of 0.04; and as Mean Absolute Percentage error (MAPE) of 3.5%. 
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galaxy-sized or cosmology-sized hydrodynamic simulations of atomic and molecular gas and 


stars where the underlying molecular cloud physicsis treated in a sub-grid fashion ( 

Pelupessy 

& Papadopoulos 

2009 

Narayanan et al. 

2011 

Lagos et al. 

2012 

Christensen et al. 

2012 ). 


Such simulations are frequently used to predict the CO-deduced versus the true H 2 gas 


distributions in galaxies (e.g. 

Pelupessy & Papadopoulos 2009 

) or make predictions of CO 

line luminosities of galaxies in a cosmological setting (e.g. 

Lagos et al. 

2012 

). Yet none 


currently takes into account the expected strong effects of CRs, driven by the often high and 
strongly evolving SFR densities found in the Early Universe. 


In the case of cosmological simulations used to predict CO SLEDs for galaxies in a 
ACDM universe (e.g. Lagos et ah 2012) knowing the temperature of CR-heated gas in the 
CO-rich and CO-poor phases is also important. In Fig|9]we show the temperatures attained 
by the molecular gas within the [% — Ccr] parameter space explored. Finally, from our 
work is now clear that Cl J=l-0, 2-1 line emission is important for recovering the CO-poor 
H 2 gas distribution along with the CO-rich one in galaxies. Cosmological simulations of 
Cl line luminosities in galaxies have been recently published (Tomassetti et al. 2014|), but 
without taking into account the important effects of CRs in regulating the CO/H 2 and Ci/H 2 
abundance ratios. 



Fig. 9.— Contours showing the corresponding gas temperatures for the parameter space 
shown in Figjl} From this Figure is also obvious that most of the temperatures attained by 
the gas are well below the temperature assumed by Bialy & Sternberg (2014) their isothermal 
treatment. 
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5. Discussion 


Quite unlike the FUV radiation fields incident on molecular clouds, CRs can destroy CO 
equally efficiently throughout their volume rather than mostly at outer cloud layers. CRs 
can thus be the main regulator of the CO-'Qisibility" of molecular gas in SF galaxies with 
the FUV-induced CO destruction a secondary one. Only for low-metallicity molecular gas 
in high-intensity radiation the FUV-induced CO destruction can surpass the CR-induced 
one. Nevertheless, given that high average FUV fields are found in places of elevated SFR 
densities, and thus also of boosted Ucr, the two effects will always operate in parallel, even as 
the average FUV radiation fields will not follow the rise of £cr to high values in a proportional 
manner in metal-rich SF environments (see 


Papadopoulos et al. 2014, for details). 


In future work we will explore the combined effects of high G 0 and Ccr cm molecular 
gas with a range of metallicities. Given the non-linear effects behind the FUV-induced 
and CR-induced CO destruction, these may be much more dramatic than either effect alone. 
FUV-intense environments with low average metallicities will be typical in the early Universe, 
either over entire galaxies (e.g. Ly-break galaxies, dwarfs) or in the outer, metal-poor, parts 
of SF disks in the process of building up their present-epoch Z(R) gradients. 


5.1. Two caveats: CR penetration, and the effects of turbulence 


Our analysis has been based on the assumption that Ucr oc Psfr, and that CRs pen¬ 
etrate into most of the volume of molecular clouds. There are good observational and 
theoretical arguments supporting both these assumptions, reviewed in Papadopoulos (2010). 
Briefly, the first assumption underlies the so-called far-IR/radio correlation established for 
galaxies whose radio and far-IR emission are powered by star formation, while the second one 
is supported by y-ray observations of H 2 -rich starbursting galactic centers, and by the ther¬ 
mal, ionization, and chemical states observed for FUV-shielded dense molecular gas cores in 
the Galaxy (Bergin & Tafalla]2007). Finally, if CRs were unable to penetrate deep inside the 
FUV-shielded regions of molecular clouds there would be no molecules present other than 
H 2 since CRs, by driving the Hjj" formation in such regions, initiate the chemical network 
beyond the basic H 2 molecule (itself formed on dust grains). 


The effects of strong turbulence on a CR-controlled CO/H 2 ratio, and its consequences 
for the CO-"visibility" of H 2 gas in SF galaxies, are difficult to assess within the static 
framework of our work. The necessity of a dynamical approach can be shown simply 
by noting that the turbulent diffusion timescales Giff are typically shorter than chemical- 
equillibrium timescales t c h, the latter approximated by the timescale needed for an Hl/H 2 
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equilibrium to be reached (e.g. [Hollenbach & Ticlens 1999). Indeed for an initially atomic 
Cold Neutral Medium (CNM) phase (and a canonical H 2 formation rate on grains of Rf — 
3 x 10 -18 (T/K) -1,/2 cm 3 s -1 ) it is t c h ~ 10' yrs, while oc l/K^m (with A'difr = (Vturb-h) ~ 
(10 24 — 10 25 )cm 2 s -1 , Xie et al. 1995) ranges from tdm ~ 3 x 10 6 yrs for quiescent GMCs 
(turbulent linewidths of Vturb ~ (1 — 5)kms -1 ) to ~ 3 x (10 4 — 10 5 ) yrs for clouds typical 


to 

of ULIRG environments (Papadopoulos et al. 2004, and references therein) The supersonic 
turbulence that dominates the velocity fields of molecular clouds is expected to influence our 
results in two distinct ways namely: a) by strongly “mixing” outer and inner cloud layers 
(and thus also the respective abundances of the various species formed in those cloud re¬ 
gions), and b) by controlling the molecular cloud mass fraction over the gas density range i.e. 
the dM gas /dn function (Vazquez-Semadeni 1994[ Padoan & Nordlund 2002, and references 
therein). Early work on the first issue by Xie et al. (1995) indicates that it is the abundances 
established on the outer cloud layers (i.e. the CO-poor/Cl-rich ones in CR-immcrsed clouds 
with realistic density profiles) that will act to “dilute” the inner, denser, CO-rich/Cl-poor 
regions. Furthermore, the typically stronger turbulence of the outer lower-density molecular 
cloud regions with respect to the inner more quiescent ones can only enhance such turbu¬ 
lent diffusion transfer of outer layer abundances inwards. These effects would then further 
reduce the CO-visibility of turbulent molecular clouds immersed in high CR energy density 
backgrounds. 


Regarding the second issue, we note that galaxies with high average Psfr values like 
ULIRGs are typically also those with much higher levels of turbulence (and hence dense gas 
mass fractions which in turn “fuel” their high Psfr)- Indeed the face-on velocity dispersion 
of the molecular gas disks in ULIRGs is er(V z )~(30-140) kms -1 (Downes & Solomon 1998), 
much higher than in isolated spiral disks (~10kms -1 ). As a result the underlying dM gas /dn 
would “shift” most of the mass of a molecular cloud towards high densities (>10 4 cm -3 ), and 
thus towards the CO-rich part of the [riu — £cr] parameter space shown in Fig|l| perhaps 
despite the also high average <(cr expected in such environments. Only hydrodynamical 
simulations of molecular clouds with such high levels of turbulence and the appropriate 
chemical networks for CO and H 2 formation and destruction can decide this issue. On 
the other hand it is worth pointing out that the cm-emission synchrotron disks of galaxies 
(rough indicators of the extend of CR propagation) are typically larger than the CO-visiblc 
molecular gas distributions. Thus starburst-originating CRs may affect the CO/H 2 ratio 
beyond the starburst region itself, where lower-density, non-SF, molecular clouds may lie. 


(Downes & Solomon 1998 
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5.2. Some observational tests 

A mainly CR-controlled CO/H 2 ratio in the molecular clouds of SF galaxies has some 
interesting observational consequences, which can be probed using the JVLA and ALMA, the 
two interferometer arrays operating in cm and mm/submm wavelengths with high sensitivi¬ 
ties. The JVLA can image the cm synchrotron radiation emitted by CR electrons (and thus 
serve as a probe of CR propagation through galactic disks and individual molecular clouds), 
while ALMA can image the spectral lines of various species emanating from CR-permeated 
molecular gas, to yield tests of a CR-controlled chemistry. More specifically: 

1. Deep imaging of the cm continuum of gas-rich SF disks along with CO J=l-0, 2-1 and 
Cl J=l-0, 2-1 line imaging can reveal the distribution of CR-permeated molecular gas, 
and probe the structural characteristics of the molecular gas distribution as recovered 
via Cl lines versus those using the two low-J CO lines. 

2. Imaging of individual molecular clouds in nearby SF galaxies in cm synchrotron con¬ 
tinuum as well as in CO J=l-0, 2-1 and Cl J=l-0, 2-1 lines in order to examine the 
effects of CRs in individual GMCs. 

3. Any molecular line observations constraining the average x(e) (ionization fraction) 
within molecular clouds, versus CR distribution (traced by cm continuum imaging), can 
reveal whether CRs penetrate throughout GMC volumes determining their chemistry 
(e.g. x(e) correlating with cm emission brightness). 

Generally in high CR energy density environments we expect H 2 clouds to be marked 
by the usual CO-bright emission emanating from the regions that are dense enough to retain 
their CO-richness, yet having Cl line emission extending much further out than those CO- 
bright regions, tracing a CO-poor phase where CO has been destroyed by CRs. Given the 
elevated temperatures of the latter phase (see Figs j7] & |9| we expect both the Cl J=l-0 
(Eio/ks 24 K) and Cl J— 2-1 (E2i/ke~62K) to be substantially excited. I 11 principle this 
allows Cl line imaging observations of SF disks to access a larger co-moving volume of the 
Universe and serve as alternative H 2 gas distribution tracers. 


6. Conclusions 

In this paper we examined the impact of CRs on the average CO/H 2 ratio in H 2 clouds. 
We find that even ordinary boosts of CR energy densities, expected in numerous SF galaxies 
in the Universe, can destroy CO (but not H 2 ), so effectively as to render potentially large H 2 
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gas reservoirs CO-poor and difficult to trace using conventional CO <7=1-0, 2-1, line imaging. 
This is a density-dependent effect, with denser regions able to retain higher CO/H 2 ratios 
than low-density ones, for a given CR energy density (and thus CR ionization rate). This 
may give rise to a second order effect with H 2 gas disks that host vigorous SF, appearing 
clumpier in CO J=l-0, 2-1 emission than they really are, their denser (and clumpier) H 2 
gas regions having retained large enough CO/H 2 ratios as to remain CO-visible, unlike the 
CO-poor lower-density ones. Furthermore, even for CR energy density backgrounds expected 
in the Milky Way, low-density gas (tt-h ~ 1 - 3 x 10 2 cm -3 ) can also be very CO-poor. The 
CR-induced destruction of CO leaves behind a Cl-rich molecular gas phase, with only small 
amounts of ClI. The Cl J=l-0, 2-1 lines are then the only remaining effective H 2 mass 
distribution tracers, encompassing both for the CO-rich and the CO-poor phase. 

Strong turbulence, often found for the molecular gas of vigorously SF galaxies, can in 
principle help maintain much of it CO-rich (and thus CO-line-visible), by placing most of 
it at high average densities. Hydrodynamical simulations of highly turbulent, CR-irradiated 
molecular gas, which incorporate the appropriate CO, H 2 formation/destruction chemical 
networks, are necessary to explore this issue. On the observational front deep Cl, low- J CO 
line, and cm continuum imaging observations of SF galaxies and individual molecular clouds 
embedded in environments of high Ucr, are key in exploring a CR-regulated CO/H 2 ratio 
in H 2 clouds. We expect the Cl line emission to extend well beyond that of low-J CO line 
emission if that is indeed so. 
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A. The numerical model 


We use 3d-pdr (Bisbas et al. 2012), a three-dimensional time-dependent astrochemistry 


code implemented at University College London and designed for treating PDRs of arbitrary 
density distribution. The code has been used in various applications examining galactic and 
extragalactic objects (i.e. Offner et al. 2013, 2014; Bisbas et al. 2014 Gaches et al. 2014). 


For the purposes of this work, we perform high resolution one-dimensional uniform-density 


clouds as discussed in Bisbas et al. (2012). We use a subset of the UMIST 2012 network 


(McElroy et al. 2013) consisting of 58 species and more than 600 reactions. The abundances 
of species relative to hydrogen we used are [He] =8.5 x 10 -2 , [C]=2.69 x 10 -4 , [O]=4.90 x 10“ 4 , 
[Mg]=3.98 x 10~ 5 , and [S]=1.32 x 10~ 5 which correspond to undepleted solar abundances 


(Asplund et al. 2009). Full technical details can be found in Bisbas et al. (2012) and in 


Bisbas et al. (2014) for further updates. 


In the simulations presented here, we use an expression for molecular line optical depths 


suitable for macroturbulent molecular clouds, 
transition ij of 


T ij — 


AijC 3 

( dv Y 

TljQi TliQj 


\d r ) VIR 

9j 


Thus in the expression of t 13 for the line 


(Al) 


(where A l3 is the Einstein coefficient, u l] is the line frequency, n*,nj the populations and g*, 
gj their corresponding statistical weights of levels i and j), we set ^ as: 


dV\ r- 

— = 0.65 

rfr /viR 


n H 


. 10 3 cm -3 


1/2 


- 1-1 


km s pc 


(A2) 


suitable for a mostly self-gravitating gas phase (Papadopoulos & Seaquist 1999), where a 


depends on the assumed density profile and % is the total H-nucleus density. This ensures 
that the optical depth of the transition ij is a local quantity everywhere within the cloud, 
as expected for macroturbulent velocity fields. For this work we adopt 0.65-y/a = 1. 
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